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Abstract. Allozyme electrophoresis was used to study the distribution of genetic variation within and among mainland 
and insular populations of the lacertid lizard Podarcis muralis from western, southern and eastern Europe. Genetic vari¬ 
ability in the species is low and genetic subdivision is high. The highest values of percent polymorphism and heterozy¬ 
gosity were found in the samples from two Tyrrhenian islands (Elba Island, La Scola Islet). The occurrence of higher 
levels of genetic variability in insular populations is probably because these populations inhabit marginal environments 
characterized by temporal-ecological instability. In these environments high heterozygosity levels can be preserved af¬ 
ter colonization events, unless founder populations are so small that bottleneck effects occur. The genetic heterogeneity 
analysis demonstrates a certain amount of genetic differentiation among local populations of P. muralis , with a relative¬ 
ly high level of genetic subdivision. Allozyme data show that genetic variation in P. muralis is distributed into two ma¬ 
jor population groups: the first includes the closely related samples from Spain and SW France, the second the geneti¬ 
cally recognizable samples from Germany, Italy, and Greece. The average genetic distance between the two groups is 
relatively high (Nei’s D = 0.059), with D ranging from 0.043 to 0.100. 

Key words. Podarcis muralis , Lacertidae, allozyme electrophoresis, population heterogeneity, Tyrrhenian islands, Eu¬ 
rope. 


INTRODUCTION 

There have been numerous surveys of the genetic struc¬ 
ture of insular populations of vertebrates, especially rep¬ 
tiles (e.g. Soule & Yang 1974; Gorman et al. 1975; Pat¬ 
ton et al. 1975). From these studies it became evident that 
many demographic, historical, and geographic factors in¬ 
fluence the pattern of genetic variation in the insular pop¬ 
ulations (e.g. Soule et al. 1973; Soule 1976). 

The Mediterranean lacertid lizards of the genus Podarcis 
seem to be particularly useful for this type of investiga¬ 
tion because they are widespread on several Mediterranean 
islands and are normally characterized by high inter- and 
intra-population morphological and genetic variability 
(e.g. Harris & Arnold 1999; Arnold & Ovenden 2002; Cor¬ 
ti & Lo Cascio 2002; Salvi et al. 2009). Although the evo¬ 
lutionary significance of the pattern of variation observed 
in these lacertid lizards has been unstudied for most taxa, 
in some cases at least it was pointed out that species which 
are characterized by a high degree of phenotypic plastic¬ 
ity in the pattern of the upper parts may have levels of ge¬ 
netic variability higher than those found in the morpho¬ 


logically low variable species (see e.g. Selander 1976; Ca¬ 
pula 1994a, 1996, 1997; Losos et al. 1997; Capula & Cec- 
carelli 2003; Caputo et al. 2008). 

In this paper, based primarily on allozyme data, the dis¬ 
tribution of genetic variation within and among mainland 
and insular populations of the lacertid lizard Podarcis mu¬ 
ralis from western, southern and eastern parts of its Eu¬ 
ropean range was estimated. Podarcis muralis was cho¬ 
sen as it is a morphologically and ecologically variable 
species occurring in a wide variety of habitats over its dis¬ 
tribution range, which extends from the northern border 
of the Iberian Peninsula to north-western Turkey, and 
throughout central and southern Europe (Arnold & 
Ovenden 2002; Corti & Lo Cascio 2002). In the northern 
part of its range this lizard is typically a thermophilous 
and lowland species, with a reduced variability in the pat¬ 
tern of the upper parts, while in the southern part it is more 
often a mountain species, occurs especially in wet and 
shady habitats, and is characterized here by high pheno¬ 
typic variability (see Capula et al. 1993, 2009; Corti et al. 
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Table 1 . Geographic and collecting data for the Podarcis muralis samples used in this study. 


Population 

Locality 

Sample size 

A 

Guadarrama (Spain) 

5 

B 

Anso (Spain) 

5 

C 

Ordesa (Spain) 

3 

D 

Deba (Spain) 

4 

E 

Albaran (Spain) 

3 

F 

Bidache (SW France) 

5 

G 

Le Chiroulet (SW France) 

5 

H 

St. Gaudens (SW France) 

5 

I 

Bonn (Germany) 

5 

J 

Cavalese (Italy) 

20 

K 

Cesena (Italy) 

20 

L 

Resceto (Italy) 

2 

M 

Chiusdino (Italy) 

6 

N 

Populonia (Italy) 

4 

O 

Uccellina Mountains (Italy) 

3 

P 

Ostia (Italy) 

10 

Q 

Elba Island, Tuscan Archipelago (Italy) 

10 

R 

Scoglietto di Portoferraio Islet, Tuscan Archipelago (Italy) 

2 

S 

Gorgona Island, Tuscan Archipelago (Italy) 

4 

T 

Pianosa Island, Tuscan Archipelago (Italy) 

2 

U 

La Scola Islet, Tuscan Archipelago (Italy) 

2 

V 

Palmaiola Island, Tuscan Archipelago (Italy) 

3 

z 

Viotia (Greece) 

6 


in press). Allozyme variation in some Italian, Spanish and 
Austrian populations of P. muralis was studied by Capu¬ 
la (1997), who provided evidence of high level of genet¬ 
ic variability in insular populations. Genetic variation and 
differentiation in the Italian populations of the species 
were recently investigated also by Caputo et al. (2008) and 
Giovannotti et al. (2010), based on the sequencing of a 
portion of a mitochondrial gene. 

MATERIAL AND METHODS 

Sampling. Samples of/) muralis used in this study were 
obtained from 17 mainland localities of western, south¬ 
ern and eastern Europe (Spain, SW France, Germany, Italy, 
Greece) and six islands of the Tuscan Archipelago in the 
Tyrrhenian Sea (Elba, Scoglietto di Portoferraio, Gorgona, 
Pianosa, La Scola, Palmaiola). The precise geographic ori¬ 
gin of each sample and the number of individuals analysed 
per population are indicated in Table 1. 


Electrophoresis. The electrophoretic analysis was under¬ 
taken for 134 specimens from all 23 localities. Standard 
horizontal starch gel electrophoresis was performed on tail 
muscle tissue, parts of which were crushed in 0.1 mL of 
distilled water. Gene products for the following 21 pre¬ 
sumptive enzyme loci were analysed: glycerol-3-phos¬ 
phate dehydrogenase (E.C. 1.1.1.8, °c Gpd), lactate dehy¬ 
drogenase (E.C. 1.1.1.27, Ldh-1, Ldh-2), malate dehydro¬ 
genase (E.C. 1.1.1.37, Mdli-1, Mdh-2 ), malic enzyme 
(E.C. 1.1.1.40, Me-1, Me-2 ), isocitrate dehydrogenase 
(E.C. 1.1.1.42, ldh-1, ldli-2), 6-phosphogluconate dehy¬ 
drogenase (E.C. 1.1.1.44, 6Pgd), glyceraldehyde-3-phos- 
phate dehydrogenase (E.C. 1.2.1.12, Gapd), superoxide 
dismutase (E.C. 1.15.1.1, Sod-1), glutamate-oxaloacetate 
transaminase (E.C. 2.6.1.1, Got-1, Got-2), creatine kinase 
(E.C. 2.13.2, Ck), adenosine deaminase (E.C. 3.5.4.4, 
Ada), carbonic anhydrase (E.C. 4.2.1.1, Co), mannose-6- 
phosphate isomerase (E.C. 5.3.1.8, Mpi), glucose-6-phos- 
phate isomerase (E.C. 5.3.1.9, Gpi), phosphoglucomutase 
(E.C. 5.4.2.2, Pgm-1, Pgm-2) (enzymes codes are accord- 
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Table 2. Chi-square values resulting from contingency y- ana¬ 
lysis of the polymorphic loci among populations of Podarcis ma- 
ralis. d.f. = degree of freedom; NS = nonsignificant. 


Locus 

No. 

of alleles 

X 2 

d.f. 

P 

Ldh-l 

2 

172.417 

22 

<0.001 

Ldh-2 

2 

33.841 

22 

<0.05 

Me-1 

4 

462.317 

66 

<0.001 

Me-2 

2 

32.622 

22 

NS 

6Pgd 

4 

222.142 

66 

<0.001 

Gapd 

2 

42.987 

22 

<0.004 

Got-1 

2 

63.179 

22 

<0.001 

Pgni-2 

2 

80.490 

22 

<0.001 

Ca 

2 

75.142 

22 

<0.001 

Gp-1 

2 

24.986 

22 

NS 

Gp-2 

2 

38.449 

22 

NS 

Gp-3 

3 

331.439 

44 

<0.001 

Total 


1580.010 

374 

<0.001 


ing to Richardson et al., 1986). In addition, three uniden¬ 
tified non-enzymatic proteins were studied: Gp-1, Gp-2, 
Gp-3. The buffer systems used, electrophoretic procedures, 
staining techniques, and loci and allele designations were 
those described by Capula (1990, 1994b). 

Analysis. Genotypic and allelic frequencies were deter¬ 
mined by direct counts from allozyme phenotypes, and the 
resulting data were analysed by various statistical meth¬ 
ods to describe the genetic structure of the P. muralis pop¬ 
ulations. Genotypic proportions expected on the basis of 
Hardy-Weinberg equilibrium were calculated by Levene’s 
formula (Levene 1949) for small samples. The statistical 
significance of departures from Hardy-Weinberg equilib¬ 
rium was estimated using a test for calculating exact sig¬ 
nificance probabilities, analogous to Fisher’s exact test 
(Elston & Forthofer 1977). To determine whether the het¬ 
erogeneity in the genotypie distribution reflects differences 
in allele frequencies, the variation in genic proportions 
among populations was subjected to a contingency y} 


Table 3. Genetic variability parameters in Podarcis muralis populations. A , mean number of alleles per locus; P, mean proporti¬ 
on of polymorphic loci; H 0 , observed mean heterozygosity; H e , expected mean heterozygosity; SE, standard error. 


Population 

A 

P 

H„ 

SE 

H e 

SE 

Guadarrama 

1,0 

0.0 

0.000 

0.000 

0.000 

0.000 

Anso 

1.0 

0.0 

0.000 

0.000 

0.000 

0.000 

Ordesa 

1.0 

4.2 

0.026 

0.026 

0.022 

0.022 

Deba 

1.0 

4.2 

0.010 

0.010 

0.010 

0.010 

Albaran 

1.0 

0.0 

0.000 

0.000 

0.000 

0.000 

Bidache 

1.1 

8.3 

0.038 

0.027 

0.038 

0.027 

Le Chiroulet 

1.0 

4.2 

0.008 

0.008 

0.008 

0.008 

St. Gaudens 

1.1 

4.2 

0.031 

0.031 

0.029 

0.029 

Bonn 

1.0 

4.2 

0.015 

0.015 

0.015 

0.015 

Cavalese 

1.0 

0.0 

0.000 

0.000 

0.000 

0.000 

Cesena 

1.1 

12.5 

0.038 

0.023 

0.040 

0.024 

Resceto 

1.1 

8.3 

0.038 

0.027 

0.042 

0.029 

Chiusdino 

1.2 

16.7 

0.064 

0.032 

0.066 

0.033 

Populonia 

1.1 

8.3 

0.038 

0.027 

0.036 

0.025 

Uccellina Mountains 

1.2 

12.5 

0.038 

0.028 

0.067 

0.039 

Ostia 

1.1 

12.5 

0.031 

0.018 

0.030 

0.017 

Elba Island 

1.1 

12.5 

0.077 

0.046 

0.055 

0.03 1 

S.to Portoferraio Islet 

1.1 

12.5 

0.019 

0.019 

0.077 

0.043 

Gorgona Island 

1.1 

8.3 

0.029 

0.021 

0.028 

0.020 

Pianosa Island 

1.1 

8.3 

0.038 

0.027 

0.042 

0.029 

La Scola Islet 

1.2 

16.7 

0.077 

0.046 

0.097 

0.046 

Palmaiola Island 

1.1 

8.3 

0.032 

0.023 

0.024 

0.017 

Viotia 

1.0 

0.0 

0.000 

0.000 

0.000 

0.000 
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Table 4. Values of Nei’s (1978) unbiased genetic distance among populations of Podarcis muralis . For geographical origin of populations (A-Z) see Table 1 . 
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analysis (Workman & Niswander 1970). The genetic vari¬ 
ability of populations was estimated using the following 
parameters: mean number of alleles per loeus (A); percent¬ 
age of polymorphic loci, at the 99% level (/’): observed 
mean heterozygosity (//„); expected mean heterozygosi¬ 
ty in Hardy-Weinberg equilibrium (H e ) (unbiased estimate; 
Nei 1978). The genetic relationships among the studied 
populations were evaluated using Nei’s unbiased genetic 
distance ( D , Nei 1978). All genetic variability and genet¬ 
ic distance measures were calculated using the computer 
program BIOSYS-2 (Swofford & Selander 1999). An es¬ 
timation of phenetic relationships among populations was 
obtained by generating a phenogram of all samples by 
means of the unweighted pair-group method with arith¬ 
metic averaging (UPGMA) based on the matrix of Nei’s 
unbiased genetic distances (Sneath & Sokal 1973). 

RESULTS 

Of the 24 electrophoretic loci analysed, ten (46%) were 
monomorphic and fixed for the same allele in all samples 
(. Mdh-1, Mdh-2, Idh-2, Gapd, Sod-!, Got-2, Mpi, Gpi, 
Pgm-1, Ada). Fourteen loci (54%) were found to be poly¬ 
morphic (oc Gpd, Ldh-l, Ldh-2, Me-l, Me-2, Idh-l, 6Pgd, 
Got-! Ck, Pgm-2, Ca, Gp-1, Gp-2, Gp-3). The Me-l and 
Pgm-2 loci only were highly polymorphic, while the oth¬ 
er 12 loci were weakly polymorphic. Four samples out of 
the 23 analysed were characterized by a unique allele (sen- 
su Slatkin 1987): St. Gaudens ( 6Pgd U0 ), Ostia (Gp-1 105 ), 
Gorgona Island (Me-2 105 ), Elba Island ( Got-l 104 ). The re¬ 
sults of the contingency x 2 analysis are given in Table 2. 
The analysis reveals that 9 out of 12 variable loci exhib¬ 
it statistically significant heterogeneity in the allele fre¬ 
quencies. This result shows that there are significant dif¬ 
ferences among the gene pools of the studied samples, in¬ 
dicating local genetic differentiation and a relatively high 
degree of substructuring among populations. Significant 
deviations from Hardy-Weinberg equilibrium in the direc¬ 
tion of heterozygote deficiences were found in the follow¬ 
ing populations and loci: Bonn (Me-l, P<0.005), Chius- 
dino (Gp-3, P<0.05). 

Genetic variability parameters (A, P, H 0 , // e ) are reported 
in Table 3. The overall number of alleles per locus (A) was 
1.07, ranging from 1.0 to 1.2. The proportion of polymor¬ 
phic loci (P) ranged from 0 (Guadarrama, Anso, Albaran, 
Cavalese, Viotia) to 16.7% (Chiusdino, La Scola Islet), av¬ 
eraging 7.25%. The observed heterozygosity (H t) ) showed 
a similar trend, ranging from 0 (Guadarrama, Anso, Al¬ 
baran, Cavalese, Viotia) to 0.077 (Elba Islands, La Scola 
Islet), and averaging 0.028. The samples from Spain, Ger¬ 
many and Greece are characterized by very low levels of 
genetic variability (Spain: average P= 1.68%, average H 0 
= 0.007; Germany: P = 3.8%, H 0 = 0.015; Greece: P = 



Fig. 1. Phenogram generated by UPGMA cluster analysis ba¬ 
sed on Nei's (1978) unbiased genetic distances among Podar¬ 
cis muralis populations. For geographic origin of populations 
(A-Z) see Table 1. 

0%, H 0 = 0) when compared with the ones from France 
(average P = 5.57%, average H 0 = 0.023) and mainland 
Italy (average / > =10.11 %, average H 0 = 0.036). Howev¬ 
er, it must be noted that Germany and Greece were rep¬ 
resented in the analysis only by a sample respectively. Per¬ 
cent polymorphism and observed heterozygosity detect¬ 
ed in island populations from the Tuscan Archipelago 
(Tyrrhenian Sea) were higher than those found in main¬ 
land samples (islands: average P = 9.52%; average H = 
0.045; mainland: average P= 5.89%; average H = 0.022); 
however the differences in polymorphism and heterozy¬ 
gosity values between mainland and insular samples were 
not statistically significant (P, P = 0.141; H 0 , P = 0.029, 
t-test). The samples from La Scola and Elba islands show 
the highest heterozygosity (H 0 = 0.077), and the sample 
from La Scola Islet exhibits the greatest genetic variation, 
with A=\.2,P= 16.7%, and H a = 0.077. 

The values of genetic distance for each pairwise compar¬ 
ison are given in Table 4. Nei’s genetic distance (D) ranges 
from 0 to 0.100, averaging 0.036. Based on the analysis 
of genetic distance data, two main population groups can 
be recognized: the first includes the samples from Spain 
and SW France, which are genetically very close (aver¬ 
age D = 0.003; D ranging from 0 to 0.010), the second 
includes all other samples (Germany, Italy, Greece) (av¬ 
erage D = 0.017; D ranging from 0 to 0.059; see Table 4). 
The ave-rage genetic distance between the two groups is 
relatively high (D = 0.059; D ranging from 0.043 to 
0.100). The comparison between the populations from 
western Europe (Spain, SW France) and Greece (Viotia) 
gives the highest genetic distances (D ranging from 0.087 
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to 0.100; see Table 4). Genetic differentiation was rather 
low among insular populations from the Tuscan Archipel¬ 
ago (average D = 0.017 , D ranging from 0 to 0.042), and 
relatively low between insular and mainland populations 
(average D = 0.040). 

The genetic relationships among the samples studied are 
presented in Figure 1. The UPGMA clustering procedure 
revealed two main clusters in the phenogram constructed 
on the basis of the matrix of Nei’s unbiased genetic dis¬ 
tances. The first cluster includes the closely related sam¬ 
ples from Spain and SW France. Within the second clus¬ 
ter the existence of four subclusters should be noted. The 
first subcluster includes the sample from Bonn (Germany), 
which is linked to the subcluster containing the closely 
grouped samples from Italy and four Tyrrhenian islands 
(Elba, Palmaiola, Pianosa, Gorgona), the third includes the 
samples from other two Tyrrhenian islands, i.e. Scogliet- 
to di Portoferraio and La Scola, the fourth contains the 
sample from Viotia (Greece). 

DISCUSSION 

The results of the allozyme analyses indicate that genet¬ 
ic variability is relatively low in P muralis. The Common 
wall lizard shows values of polymorphism and heterozy¬ 
gosity higher than those estimated by Capula (2004) for 
P. raffonei ( P = 4.8%; H= 0.011), similar to those ob¬ 
served by Capula & Ceccarelli (2003) for Italian popula¬ 
tions of P. sicula (P = 10%; H 0 = 0.029), and lower than 
(i) those detected in the phylogenetically related P wa- 
gleriana from Sicily (P = 15%; H 0 = 0.037; Capula, 1994b) 
and P tiliguerta from Sardinia and Corsica (P= 22%; H ( = 
0.066; Capula, 1996), (ii) the average ones calculated by 
Capula (1990) for nine species of the genus Podarcis (P 
= 13%; H 0 — 0.053), and (iii) the average ones calculated 
by Nevo (1978) for 17 species of reptiles (P = 22%; H a 
= 0.047). The highest values of heterozygosity were found 
in the samples from Elba Island and La Scola Islet (Tus¬ 
can Archipleago, Tyrrhenian Sea), whereas the lowest ones 
were observed in some samples from Spain (Guadarrama, 
Anso, Albaran), Italy (Cavalese) and Greece (Viotia). 
Based on the theory (see e.g. Nei et al. 1975, Gorman et 
al. 1975, 1978) we expected to find low levels of genet¬ 
ic variability in the insular samples of P. muralis. Our re¬ 
sults were in some way not congruent with these expec¬ 
tations, as some samples from the Tuscan Archipelago 
(e.g. Elba, La Scola) were characterized by levels of per¬ 
cent polymorphism and heterozygosity higher than those 
found in most of the populations from mainland Italy, 
Spain, SW France, Germany and Greece. This result in 
agreement with the allozyme data provided by Capula 
(1997) indicating that the insular P. muralis populations 
from Elba Island, Isolotto di Porto Ercole Islet and Argen- 


tarola Islet (Tyrrhenian islands) are characterized by lev¬ 
els of percent polymorphism and heterozygosity higher 
than those found in the populations from the Italian Penin¬ 
sula, and much higher than those observed in the conti¬ 
nental populations from the Italian Alps, Spain and Aus¬ 
tria (Capula 1997). Within the lacertid lizards, levels of 
genetic variation are known to be large only in mainland 
populations of a few species (e.g. Adriatic populations of 
P. sicula : average H = 0.09 (Gorman et al. 1975); Acan- 
thodactylns spp.: average H= 0.18, average P = 50%; 
Blanc & Cariou 1980), while populations living on relict 
islands and on tiny fringing islands (i.e. very small islands 
that are separated by a short linear distance from the moth¬ 
er island or continent) are usually characterized by very 
low values of per cent polymorphism and heterozygosi¬ 
ty (Gorman et al. 1975). However, the investigated 
Tyrrhenian islands can be considered as relict islands, as 
their lizard populations are genetically differentiated from 
the mainland ones. Moreover, one of the islands consid¬ 
ered here (Elba) is a large island (223 km 2 ), while the oth¬ 
er (La Scola) is a tiny island (0.016 km 2 ), and both are 
separated by a relatively short geographic distance from 
mainland (peninsular Italy). 

Among reptiles, high levels of genetic variation found in 
populations of some species (e.g. Podarcis sicula, Cne- 
midophorus tigris) arc ascribed to high vagility and con¬ 
sequent low levels of inbreeding (Gorman et al. 1977). 
This does not seem to be true in the genetically highly vari¬ 
able species of Acanthodactylus, as these are territorial 
lizards (Blanc & Cariou 1980). Podarcis muralis is a ter¬ 
ritorial lizard as well (Steward 1965), but in this case on¬ 
ly some populations — almost exclusively island popula¬ 
tions (Elba and La Scola: this paper; Elba, Isolotto di Por¬ 
to Ercole, Argentarola: Capula 1997) - exhibit high ge¬ 
netic variability. As suggested by Capula (1997), this is 
probably because the insular populations inhabit margin¬ 
al environments characterized by temporal and ecologi¬ 
cal instability. According to Lewontin (1974), in such en¬ 
vironments no particular genotype is favoured for long pe¬ 
riods and natural populations usually show levels of ge¬ 
netic variability higher than those found in more stable en¬ 
vironments. On the basis of these considerations, finding 
greater genetic variability in insular populations of P. mn- 
ralis could indicate that high heterozygosity levels can be 
preserved after colonization events in marginal popula¬ 
tions of vertebrates, unless founder populations are so 
small that bottleneck effects occur. 

The genetic heterogeneity analysis demonstrates a certain 
amount of genetic differentiation among local populations 
of P. muralis, with a relatively high level of genetic sub¬ 
division. Allozyme data show that, at the scale of the study, 
genetic variation in P. muralis is distributed into two ma¬ 
jor population groups: the first includes the closely relat- 
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ed samples from Spain and SW France, the second the ge¬ 
netically recognizable samples from Germany, Italy, and 
Greece. Genetic distance values found between the two 
groups are relatively high (Nei’s D ranging from 0.043 to 
0.100), although falling below those normally encountered 
comparing populations of well recognized biological 
species of the genus Podarcis (see e.g. Mayer 1981; Thor¬ 
pe 1983; Capula 1994a, b, c, 1996). The high genetic affin¬ 
ity between the French and Spanish samples is congruent 
with their geographic origin, as French samples are from 
Pyrenean localities (F-H) close to the Spanish ones (A-E). 
On the other hand, the large geographic distances sepa¬ 
rating the French localities from the German one (Bonn) 
could explain the relatively high genetic differentiation ob¬ 
served between the samples from these countries, which 
cluster separately in the UPGMA phenogram (see Fig. I). 

The data presented here are in agreement with the results 
of the allozyme investigations earned out by Capula 
(1997) on some P. muralis populations from Italy, Spain 
and Austria, and are congruent with the residts of the mo¬ 
lecular investigations (analysis of mitochondrial DNA se¬ 
quences) carried out by Caputo et al. (2008) and Giovan- 
notti et al. (2010) on several Italian samples, which indi¬ 
cate a certain amount of molecular divergence among P. 
muralis populations, and a pronounced geographical struc¬ 
ture of the Italian populations. 
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